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Abstract 

We describe a mathematical model of digestion in the small intestine. The main interest 
of our work is to consider simultaneously the different aspects of digestion i.e. transport 
of the bolus all along the intestine, feedstuffs degradation according to the enzymes 
and local physical conditions, and nutrients absorption. A system of coupled ordinary 
differential equations is used to model these phenomena. The major unknowns of this 
system are the position of the bolus and its composition. This system of equations is 
solved numerically. We present several numerical computations for the degradation, 
absorption and transport of the bolus with acceptable accuracy regarding the overall 
behavior of the model and also when challenged versus experimental data. The main 
feature and interest of this model are its genericity. Even if we are at an early stage of 
development, our approach can be adapted to deal with contrasted feedstuffs in non- 
ruminant animal to predict the composition and velocity of bolus in the small intestine. 

Keywords: Digestion, Small Intestine, Modeling, Ordinary Differential Equations, 
Enzymatic Degradation, Transport 



1. Introduction 

The main step of digestion and absorption along the gastrointestinal tract takes 
place in the small intestine for non ruminant animals. To reach an optimized composi- 
tion of available nutrients due to their behavior in the digestive system, it is necessary 
to understand and predict the digestion and absorption of the ingested feedstuffs in the 
small intestine |4, 8 6|. It is also now well-known that the use of implanted experi- 
mental devices may modify the dynamic of digestion in the small intestine ifTOllSTI . 

Several models have been developed representing the digestion and transport of 
bolus in the small intestine. In the model of [IJ digestion and absorption are integrated 
and represented considering only the polymers and individual absorbable end prod- 
ucts. The transit through the small intestine is modeled mainly as a result of gastric 



Preprint submitted to Theoretical biology 



Octobers, 2011 



emptying. No peristaltic wave is taken into account, and the bolus contained only the 
dry matter. |5| describe the digestion and absorption using the plug flow reactors to 
encapsulate complex digestion phenomena in a simple set of equations. Different rate 
of absorption and degradation are involved: first order kinetics, Michaelis-Menten and 
the sigmoid ones. A detailed model of the intestinal propulsion is provided by Il9l l20l . 
However, these models portray the transport of bolus simplistically, or they represent 
only a limited number of different processes involved in digestion. 

This article tries to go further in the modeling of digestion in the small intestine 
by considering the different steps of digestion i.e. the transport of the bolus all along 
the intestine, feedstuffs degradation according to the enzymes and local physical con- 
ditions, and nutrient absorption. Therefore a system of coupled ordinary differential 
equations is used. The major unknowns of this system are the position of the bolus and 
its composition. 

In fact, several models are presented reflecting the modeling process at its different 
stages with our attempts to make it more realistic by inclusion of more sophisticated 
and relevant biological phenomena and chemical transformations. We decided to de- 
scribe the different steps with the assumptions leading us to our choices instead of 
presenting only the last model since the whole process by itself may help to under- 
line relevant questions to be further discussed. Of course, this modeling process is an 
iterative one and is still going ahead in directions which are described in Section[6] 

Our models intend to be a mechanistic approach of feedstuffs digestion even though 
simplifications were included according to participatory approaches between biologists 
and mathematicians. However they involve a lot of different unknowns and parame- 
ters, and require a numerical software to obtain suitable approximation of the solutions 
since it is hopeless to obtain explicit ones. Scilab software was used to perform these 
numerical computations [^. 

In all our models, we try to estimate the parameters using data from scientific lit- 
erature. When these data are not available, we assume the reasonable values for the 
parameters. 

The article is organized as follows : Section 2 is devoted to present the main as- 
sumptions of our models and most of our notations. In Section 3, we describe the 
transport equations. In our four different models, we point out that there are only two 
different ways of modeling the transport of the bolus in the intestinal tract. The main 
differences concern the degradation itself, with several possibilities for the composi- 
tion of the bolus, for the enzymatic reactions and the water influence. The outcoming 
stages (4 different models) are presented in Section 4, with the key assumptions and 
characteristics of each model. Section 5 is a comparison of these four models and of the 
numerical results of the most sophisticated model (Model 3) versus some experimental 
data from the literature. Finally, in Section 6, we criticize our models and describe the 
perspective. 



The reader can perform its own numerical experiments, witli tlie possibility of changing the parameters, 
using our Scilab softwai'e online at the URL : http://www.lmpt.univ-tOurs.fr/modingre 
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2. General Hypothesis and Synthetic Presentation of the Different Models 

Common assumptions to all models are the following ones. 

(i) The first simplification concerns the small intestine representation itself In- 
stead of taking into account its complex geometry, it is represented as a one- 
dimensional interval [0, L\. The position of the bolus in the small intestine at 
time t is given by x{t) e [0, L] (cf. figur^. The equation of the transport of 
bolus along the small intestine is described on x{t). 

{ii) The bolus is treated as a homogeneous cylinder with a fixed length i and variable 
radius R{t) which is a function of time t. To locate this cylinder, we use the 
position x{t) of its center. This assumption is justified by the general shape of the 
small intestine's segments as well as the observation of the real bolus in animals' 
small intestines. As the length of bolus is assumed fixed, the degradation only 
changes the radius. Degradation of substrates is obtained by enzymatic reactions 
with possible subsequent absorption by the intestinal wall 1 10 1. 

{Hi) The enzymes which participate in enzymatic hydrolysis included in our models 
consist systematically in pancreatic and brush-border ones with the possibility of 
including exogenous and gastric ones. The enzymes' activity depends on the pH 
of the small intestine at each point along its length. The brush-border enzymes on 
the intestinal wall are assumed to be always in excess ifTOl . 




Figure 1: The Cylindric bolus with its different characteristics 

The bolus moves through the intestinal tract because of the pulses resulting from 
peristaltic waves and gastric emptying which gives an initial velocity to the bolus I.IOJ . 
Peristalsis are series of wavelike contractions occurring in the smooth muscle layer of 
the gastrointestinal tract. It is a physiological process that results in intestinal motility 
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and propulsion of ingested food along the intestine. It starts as a ring-like constriction 
initially which later moves mostly forward along the intestine. Moreover it might be 
assumed that it helps the bolus to be digested by spreading the food particles along the 
intestinal wall for effective digestion and absorption. Therefore to model the transport, 
bolus movement is connected to pulses all along the intestine with the initial velocity 
coming from the gastric emptying effect. Efficiency of pulses is proportional to the 
volume of the bolus and it is inversely correlated to the distance between the bolus 
and the pylorus ||3]|T8l. The bolus speed is assumed to be slowed down by the local 
conditions in the small intestine lumen (friction on the borders, viscosity effects, ...etc). 
The effects of these different local factors depends on the composition of the bolus, 
and in particular its dilution. 

For the bolus content, the following assumptions and notations are used throughout 
this paper 

• The bolus includes a single species whose total mass is denoted by A. In the most 
sophisticated model presented in this work, A is composed of As, Ans and And, 
in other words A =As+Ans + A^^ . A^ is the mass of solubilized fraction of A 
which can be hydrolyzed in the presence of the enzymes. The index s stands for 
"solubilized". Ans is the mass of "non-solubilized" fraction of A, for example 
dry starch or the non-emulsified lipids. Transformation of Ans into As requires 
a sufficient quantity of water. Regarding lipids emulsion, we assume that the 
bile salts are in excess. The mechanism Ans ^ is described through an 
equilibrium property depending on the water quantity in the bolus. Finally And 
is the mass of non-degradable A, which enters and leaves the small intestine 
without any change. For example the vegetal fiber in feed are poorly digestible. 
Moreover the fiber matrix of feedstuffs or the anti-nutritional factor content can 
be responsible for a reduction in the digestibility of some amino acids in some 
feedstuffs |17|. 

• The quantity B is the mass of product obtained from As by enzymatic reactions, 
it is composed of Bint and Babs 



The quantity B^nt is the product of hydrolysis due to gastric and pancreatic en- 
zymes, the index int stands for "intermediate" substrate which is not yet ab- 
sorbable since it is not fully degraded. This transformation has to be completed 
by a second one at the border of the small intestine via the brush-border enzymes 
(e.g. : proteins being degraded to polypeptides and afterwards to dipeptide or 
amino acids, which are absorbable). This second transformation is also able to 
give Babs directly from A. The quantity Babs is the absorbable fraction with 
index abs indicating "absorbable". 

• The quantity e represents the gastric and pancreatic enzymes. 

• The quantity W is the mass of water in the bolus and [W] indicates the proportion 
of water in the bolus : W/{A + B + W). 
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• The quantity V{t) denotes the volume of bolus which is equal to (A+B + W) /p, 
where p denotes the density. For the sake of simplicity, we assume that all the 
substrates of bolus have the same density p. The total mass of the bolus is {A + 
B + W) (t) at each moment. 

Digestion consists in the transformation of digesta to absorbable nutrients through en- 
zymatic hydrolysis. Volumic transformation is the degradation of Ag into Bi^t inside 
the bolus and transformation on the bolus surface is the degradation of both Ag and 
Bint into Babs on a thin layer around the bolus. The following hypothesis are added 
progressively with upgraded versions of the model. 

Model 1. In the first model, the whole bolus is considered to be solubilized (A =As). A is 
hydrolyzed thanks to gastric enzymes and becomes directly absorbable nutrients 
(B =Babs )■ In this model, brush-border or pancreatic ones are not taken into 
account. Such mechanisms are associated for example with the consumption of 
disaccharides (resp. monosaccharides) such as milk sugar(resp. glucose). 

Model 2. This model is an attempt to have a more realistic modeling of degradation. The 
bolus is still assumed to be completely solubilized. The absorbable nutrients can 
be obtained by two ways: either by a direct surfacic transformation A Bats or 
through a first volumic degradation A Bin followed by a second one Bint — > 
Babs at the bolus surface by brush-border enzymes. 

Model 3. This model includes the solubilization of the bolus in presence of water. A is 
splitted into Ag, Ans , And ■ Equations are added to express the equilibrium 
As ^ Ans which depends on the quantity of water. The non-degradable part of 
bolus enters and leaves the small intestine without any mechanical or chemical 
change in its initial form. A key feature of this model concerns the transport 
of bolus along the small intestine since it is connected to the quantity of water 
through lubrification effects. 

Model 4. This model is a simplification of the previous one by mathematical arguments. 

Through homogenization methods it is shown that the acceleration can be aver- 
aged and an equation with this averaged acceleration is substituted for the pulses 
in transport equation. Detailed models are described in Section 4. 



3. Transport 

We present a mathematical formulation of the transport of bolus in the small in- 
testine. It is based on the physiology of the pig's small intestine to get consistent pa- 
rameters. The duodenum is characterized by oscillatory electrical events (slow waves) 
occurring at a rate of 18/min. After food intake, some of these waves are associated 
with spikes bursts which are responsible for contractions and therefore propelling the 
bolus through the small intestine. We assume only 6 of the 18 slow waves by minute 
are followed by the spikes. It leads to one efficient contraction every 10 seconds ID. 
The mean transit time of each peristaltic wave is assumed to be 150 minutes to move 
along small intestine from duodenum to the end of ileum according to |7 1 and 1 14] . We 
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assume also that the pig's small intestine is about 18 meters (12], hence the average 
velocity c of these waves is 7, 2 m/h. Each peristaltic wave takes x{t)/c seconds to 
reach the bolus in position x{t), therefore the pulse which pushes the bolus in time t is 
generated in duodenum at time t — x{t) jc. 

If v{t) denotes the velocity of the bolus (v{t) = "j^(^))' ^^he effect of pulses is 
modeled through the following equation 

d 

the term — y represents the pulses which are defined as a periodic function of period 

at 

10 seconds such that 

"1° d 

and for t < 0, we assume y {t) ~ 0. 

Over a period, each pulse is an approximation of a Dirac mass of the origin. There- 
fore we define it as a function with the value 1/e during a very short interval of time e 
and at all other time. 

According to ifTSl and ifTSi the efficiency of the peristaltic waves increases with 
the size of the bolus and decreases with the distance from pylorus. We assume that all 
these dependences are affine, namely 

d^x , . d , , , ^ , M Co + ciVit) 

where, cq and ci are determined under the assumption that the acceleration depends 
linearly on V{t). THe constants a and b are obtained from experimental data. 

The intestinal lumen is a confined environment which prevents the bolus to move 
perfectly according to the previous equation : the bolus has to work its way through the 
small intestine and is also submitted to the friction with the intestinal wall. All these 
friction effects are related to the "viscosity" of the bolus and we have two different 
ways to model the friction term : either as a constant effect which is independent of 
the bolus composition (models 1 & 2) or with a lubrification effect coming from the 
proportion of water in the bolus (models 3 & 4). More specifically, in Equation ([T]l 
below, the coefficient K{t) is either constant in models 1 & 2 or K/[W] in models 3 
& 4, where X is a constant. 

The final equation of transport is therefore the following one 

d^x , , d , , , X , M cn + CiV(t) .dx,. 
^(0) =vo, x{0) = 

4. Digestion 

Digestion is a mechanical and chemical process by which the feedstuffs molecules 
are broken down to the smaller ones by enzymes in order to to make them available 
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for absorption. The uptake of the obtained nutrients is mjiinly by absorption. In this 
section the different steps of modehng are detailed. 



4.1. Model 1 

In this model the bolus is assumed to be completely solubiUzed (^=^4^). We also 
assume that the necessary enzymes for hydrolysis are mixed with the bolus in the stom- 
ach. The product of following reaction is directiy absorbable {B =Babs ) 

A + e^ B. 

e denotes the gastric enzymes. The first aim of this model is to define the variation 
of the bolus which means the amount of the different substrates A, B and e at every 
time in the cylinder, and the second is to locate the bolus along the small intestine. We 
assume that, the evolution of A or its volumic transformation depends on its mass at 
each moment and the enzyme activity. This equation foUows the law of mass action 

^ = -Ck(x,e)A 
at 

where C denotes the degradation rate and, k{x,e) is the enzyme activity which depends 
on the pH of the small intestine and the presence of the enzymes at each point along it. 

The product B of the volumic transformation of A is absorbed by intestinal wall 
with a constant rate kabs 

dB 

— = C'k{x, e)A - kabsB. 

There are also the degradation and inactivation of the enzymes along the small intestine 

de 
dt 

where kg is the rate of degradation of the enzymes which depends on their types. The 
activity of each enzyme as a function of pH of small intestine is roughly known. We 
know also the pH of each point along the small intestine. The composition of these 
two functions gives the enzyme activity at each point x along it. 



4.2. Model! 

In this second model, the presence of pancreatic enzymes in the small intestine as 
well as the brush-border ones on its wall are considered. The pancreatic secretions help 
neutralizing the stomach acid as they enter the small intestine. They also contain pan- 
creatic enzymes. The level of the secretions is a function of volume and composition 
of the bolus entered the small intestine. The brush border enzymes are the enzymes 
for the terminal stage of digestion which is the surfacic hydrolysis. Contrary to the 
pancreatic enzymes they are not free in the intestinal lumen, but rather, in the plasma 
membrane of the enterocyte. 

We assume that the bolus is completely solubilized. The product of the hydrolysis 
B consists in B^ and Babs (B =Bint +Babs )• 
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The following scheme represents the chemical reactions of the bolus in this model 

Binti Babsi Bint B^bs- 

The first reaction takes place inside the bolus by pancreatic and gastric enzymes, the 
second and the third ones take place on the surface of the bolus. 

The degradation of A in this model is the result of the volumic hydrolysis of A as 
in Model 1, and its surfacic hydrolysis by brush-border enzymes 



dA A 

— = -Ckix,e)A-Cabs{27rR£) 

at A + Hint + i'abs 

= -Ck{x,e)A-2CabsVT^P ^ 



{A + Bint + B^bsy/^' 

the second term represents surfacic transformation of A to Babs ■ We recall that the 
mass of the bolus in this model is 

A{t) + Bintit) + Babsit) = pV{t) = fmR^{t)l, 

and therefore the lateral surface of the cylinder is given by 

27Tm = 2y/7Tl/p{A + Bint + BabsY'"^ ■ 

This transformation depends on the fraction of A on the surface of the bolus which 
is written by {2TrR£) — — . The unit of the degradation coefficient per unit 

A + Bint + Babs 

of surface and time, Cabs , is g.m~'^ .s~^ . 

After a distance traveled by bolus of about 5% of the total length of the small 
intestine which is approximatively 85 cm in an growing pigs, the input of secretions 
starts and it stops after a distance of a meters traveled by bolus. We assume their mass 
is about /3% of the bolus mass. In the following equation, the effect of these secretions 
on the variation of A is taken into account 

rl A 1 rlr 

— = ... + ln(l./3)--x((x(s)-0.85)/a)A 

where x is a localization function in the above equation which reflects the fact that se- 
cretions arrive in the small segment of the intestine, say between 0.85 cm and 0.85 + a 
cm. 

The product of volumic hydrolysis. Bint > participates in the creation of Babs on 
the surface of the bolus. Therefore its variation is modeled by 



^ = Ckix,e)A + ln{1.25)^^xiixis)-0M)/a)Bint 

B,,u 



-2Ciabs\/'^i/P 



{A + Bint + Babsy/^' 
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The absorbable nutrients on the bolus are not absorbed instantaneously (5). In this 
model we assume that the absorption rate follows Michaelis-Menten mechanism. The 
constant kabs is the maximal rate of absorption at saturation, k is the Michaelis constant 
which is half saturation 

2^/^ 



dt ^ "^{A + B.nt+BabsV/^ ""'k + Babs' 

4.3. Models 

In this model the ingested food consists in Ans, And, As and water {A — 
And + + W). We incorporate two effects of water on digestion : the first one is the 
dilution of the bolus and its impacts on degradation and absorption and the second one 
is the lubrification and its consequences on the transport. 

We assume that the evolution of Ag and Ans aims at reaching an equilibrium in 
which the ratio between Ag and Ans is fixed and depends only on the proportion of 
water, namely As = iJ.{[W])Ans stressing that solubilization of Ans depends on 
bolus dilution. From the mathematical standpoint, we write this evolution as 

^=-ks(^li{[W])Ans-As^, (2) 

where /i is a linear function of water and the constant kg represents the return rate to 
equilibrium. 

The amount of water in the intestinal lumen is regulated by several complex bio- 
logical phenomena. In fact the proportion of water in the bolus aims at reaching [Wo] 
in a rather fast way which we translate it on a mathematical standpoint 

^ = -kUm ~ [W,]) + ln(l./3)-^x((a:(5) - 0.85)/a)[W] (3) 
at a at 

where k^ is large enough to reach the equilibrium in an adequate time. The second 
term of above equation is the fraction of water in pancreatic secretions. 

The variation of As depends on its degradation by volumic and surfacic hydrolysis, 
and contribution of pancreatic secretions as in previous models. It also depends on the 
equilibrium with Ans resulting from equation Q which is the first term of equation 
below 

^ = ks(^^i {[W]) Ans - As^ - Ck{x, e)As{t) 

- 2C.bsVWp^^^ ^ ^ [W] (4) 

1 dr 

+ ln(1.25)--x((2:W-0.85)/a)A,. 
a dt 

The variation of absorbable nutrients depends on the creation of Babs by enzymatic 
hydrolysis of As and Bint and its absorption by intestinal wall 
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I T"/ ^abs-^ ~(~ Ciabs-^int v-tTri 



dt ^ {As+ Ans + And + B,nt +W + BabsV^^ 

, Babs 

k + Babs 

The non-degradable fraction of A, namely And, enters in the small intestine and leaves 
it without any change in its structure. 

As we already indicated in Section [3] lubrification of the bolus depends on the 
presence of water. For this model, the friction coefficient in equation (|6]l is written as 

Kit) = T^^. 

4.4. Model 4 

This model is a mathematical simplification of the transport equation by means 
of homogenization theory. Homogenization theory is concerned with equations with 
rapidly oscillating coefficients and its aim is to provide an "homogenized" or "aver- 
aged" equation which is a limiting equation when the frequency of the oscillations 
tends to infinity. Advantages of homogenization theory are clear : on most occasions, 
it is simpler to use homogenized equations (for example to compute the solution) and, 
when the frequency of oscillations is over a given value, this approximation of the real 
equation by the homogenized one may be rather accurate as seen in the next section. 

Homogenization problems for ODEs were studied by 1 11 1 but it is worth pointing 
out that our particular case does not fall into the theory described in 1 11 1. Fortunately 
the specific structure of the transport equation allows us to do a complete analysis of 
the problem and even to compute explicitly the averaged equation. 

More specifically, in the transport equation, pulses reach the bolus every 10 seconds 
approximately. Compared to the time scale of digestion phenomena (the bolus stays 
in the small intestine for several hours 1211 this represents a very high frequency 
and causes very rapid variations in the velocity of the bolus (see the velocity profile in 
Figure [3]l. 

We can prove mathematically that the pulses can be averaged out in an appropriate 
way and we can replace the rapidly varying velocity by a slowly one. 

In the simplest case, by normalizing the pulses, we assume that their mean effect 
over a period is e(e). Thus, over a time t = Ne, their mean effect is Ne{e) = te(e) /e. 
On the other hand we assume 

lim e(e)/e = r 
the mean effect over a time t is therefore 

lim A^efe) ie(e)/e = tr. 

Inserting this equality in transport equation ([T]l, the homogenized transport equation 
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reads 



dx 



(0) 



it) 



vo, x{0)^0 



dit) 



co+CiV{t) 
a + bx(t) 



[W]{t) dt 



K dx 



it) 



where a is the averaged effect of the pulses 



Its value is 



d[t) t(1 



1 dx 
c dt 



5. Results 

In the first part of this section, the graph of degradation of model 4, and the graph of 
transport of model 3 and 4 are developed. The second part concerns the evaluation of 
the last model by comparing its outputs with experimental data. Only a hmited num- 
ber of outputs can be compared because of the lack of experimental data. However, 
the model is evaluated in relation to our objective which is developing a mathemati- 
cal model that takes into account the physiology of the small intestine and process of 
digestion in it. 

5.1. Digestion 

The graph of digestion of model 4 is shown in figure |2] We should at first initiate 
the bolus composition. These initial conditions vary following the different types of 
feedstuffs. We fixed the initial value of Ans as three times that of Ag. We dilute Ans 
by two times its volume water Solubilized substrate A^ and non-solubilized one Ans 
reach a dynamic equilibrium all along the small intestine, as explained in Section |4] 
This balance is reached rapidly at the beginning of the small intestine due to the large 
difference in quantity between these two substrates. The result of this equilibrium is 
the increases of the value of As and the decreases of the values of Ans, as seen in the 
graph of digestion. The inverse process might take place by lack of water. The absorp- 
tion curve corresponds to the collected absorbable nutrients from x = to a; = x{t), 
where x indicates the location of bolus in the small intestine. Obviously, the graph of 
the fourth model contains more details about different steps of digestion than the first 
two graphs thanks to the model structures. The digestion graph of model 4 is similar 
that of model 3. The only change in model 4 deals with the transport equation. 
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Figure 2: Digestion through Model 4 



5.2. Velocity 

Figure |3] provides a numerical evidence of the homogenization phenomena. The 
graph of transport resulting from models 3 and 4 is shown in figure |3] The effect of 
pulses on the curve of velocity is obvious. However, using the homogenization theory 
in Model 4, we obtain a smooth graph of velocity which replaced that of Model 3. 
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Figure 3: Velocity of the bolus versus Time 



5.3. Model Evaluation 

For a specific family of nutrients, here starch, digestion is calculated using Model 
4 and is compared to data reported by To parameterize adequately the model, 
we adapt the enzyme activity of the last model to the activity of amylase in the small 
intestine. Amylase is the enzyme required for degradation of starch. The optimal 
activity of pancreatic amylase is in neutral pH [TOl. 

Pancreatic secretions have no impact on the variation of Ag since there is no starch 
from this source. 

The inputs of model are only Ans and W which are Starch and Water. The outputs 
are the values of these substrates at the end of ileum. The data in the article of [2 1 are 
for purified protein free wheat starch, agreeing with our hypothesis for the composition 
of the bolus [T| The outputs concerns the collected data after at the end of ileum. 
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Table 1 : Digestion of Starch in Pigs by Modeling : Comparison between simulated and 
experimental data by 121 





Experimentation 


Modeling 




Input(g) 


Output(%) 


Input(g) 


Output(%) 


wet digesta 


2571 


8 


113.10 


5.33 


dry matter 


688 


0.50(g) 


37.70 


0.04 



Regarding the data presented table [T] percentages of dry matter and wet digesta 
collected at the end of the ileum are approximatively the same as the output of model 
4. The difference between inputs is due to the simulation calibration which takes into 
account only one bolus i.e. a fraction of the daily meal. However, differences between 
outputs are low in percentage enabling to conclude that the model can roughly simulate 
very simple situations. 

5.4. Sensitivity Analysis 

Sensitivity analysis is performed to identify the key parameters affecting the diges- 
tion process. The chosen parameters are set at 5% and 50% of their original values. 



Output 


Parameters 


As 


C, Cabs 


Babs 


Cabsy Ciabsi ^abs 


V 


a,b, Co, Ci, K 



Studied digestion parameters are C and Cabs for degradation of Ag, and Cabs, Ciabs 
and kabs for the absorption of Babs- 

If y is the output and 6 the parameter, the relative variation of y can be expressed 
as follows 

\ye - ye+Ae\ 

ye 

5.4.1. Influence on Ag 

Both parameters C and Cabs are overestimated by 5 and 50%. The figure |4] shows 
the relative variation of Ag at each moment. The relative variation of Ag resulting from 
5 and 50% values of C is not meaningful. The parameter Cabs has the largest effect on 
Ag degradation. Observing the graph of relative variation of Ag, figure|4] we conclude 
that increasing the value of Cabs increases the relative variation value with time. . 
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Figure 4: Relative variation of As regarding to C, Cabs 



5.4.2. Influence on Babs 

The parameters kt. Cabs and Ciabs are overestimated by 5 and 50%. The figure |5] 
shows the relative variation of Babs by time. The quantity of Bint being very small in 
the model, the effect of changing the parameter Ciabs is neglectible on the relative vari- 
ation of Babs by time. The quantity Babs is very sensitive to the variation of parameter 
Cabs firstly because of the high quantity of Ag, then its influence decreases because of 
decreasing quantity of Ag over time. The quantity Babs is dependent on kt because of 
the large impact of kt on the nutrient absorption rate. 
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Figure 5: Relative variation of Babs regarding to Cabs, Ciabs and kt 



6. Conclusion and Perspectives 

This model is obtained from simplified biological assumptions and it can be used to 
illustrate generically the rate of degradation and absorption all along the small intestine. 
This is a global model of digestion of a bolus composed of one substrate and water. This 
section is devoted to a discussion on the current state of our modeling, our assumptions, 
difficulties and on the future development of the model. 

The first assumption to be discussed is the "cyUnder" one. As mentioned earlier, 
it was introduced for technical reasons. Our aim was to solve a partial differential 
equation with very different scales of times (pulses arising every 10 seconds while the 
whole digestion process in the small intestine lasts for several hours 11211 121') and with an 
highly variable domain with contrasted scales (few centimeters for the bolus compared 
to the 18 meters of the small intestine). Solving this PDE seems unreasonable since it 
was leading to the usual diffusive phenomena and large errors. We also notice that these 
18 meters of the small intestine were rather empty and therefore we were spending a 
lot of time to compute functions which were very often 0. 
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From the transport standpoint, the "cyHnder" assumption can be seen as a La- 
grangian method, the ordinary differential equations on x{-) being (essentially) the 
characteristic curves of the transport equation. This is the first justification of this hy- 
pothesis, the second being the direct observation of animals bolus which convinces us 
that it can be represented as a cylinder, even if its geometrical characteristics could be 
more complicated. However we have to work more on the evolution of the length of 
the cylinder 

A more Eulerian "compartmental approach" is studied simultaneously but we are 
still facing difficulties for modeling the transport and, specifically to reproduce some 

features of the "cylinder" model. 

The transport equation seems to take into account rather closely the phenomena 

Co + ciA 

which are described by the experts . It will be difficult to validate the term ; — 

a + bx 

and to have a precise idea of the value of the different constants but such a modeling 
seems more appropriate than trying to use a complicated fluid mechanics approach 
whose laws may not be valid in this very confined domain. The same remarks hold for 
the effects of the water : it seems correct even if a relevant validation will be difficult. 

For food digestion and absorption, we are only at a first stage of modeling. The 
absorption phenomena were not studied expUcitly leading to required further develop- 
ment with a focus on the assumed interactions between the animal physiological status 
and absorption. The spatial aspects (location of the absorption) were clearly neglected 
so far. 

For digestion, the next step will be to mix different nutrients and adapt the enzyme 
breakdown to each of them. We have also to examine more closely the respective 
effects of the different categories of enzymes together with the role of the water. More- 
over interactions between nutrients on the digestion processes should be questioned. 

As a conclusion of this first stage of modeling, consistent behaviors of the model 
were reached. Moreover, the simphcity of the current model allows easy developments 
in any directions. Our next target will therefore be to iterate the model development 
according to the above proposed research areas. 
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